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Abstract 


This study is concerned with the structure of pulsar magneto- 
spheres and the acceleration mechanism for charged particles in the 
magnetosphere. We follow the pulsar model developed by P. A. Sturrock 
(I 97 I) and asstime that charged particles are accelerated from each 
polar cap of a pulsar. These particles produce gamma rays via curvature 
radiation which in turn produce electron-positron pairs which are 
ultimately responsible for the observed radio emission. This model 
requires large acceleration of the particles near the surface of the 
star. 

The required acceleration has not been produced in earlier pulsar 
models. We have developed a theorem which shows that particle accelera- 
tion cannot be expected when the angle between the magnetic field lines 
and the rotation axis is constant (e.g. radial field lines). If this 
angle is not constant, however, acceleration must occur. 

We have investigated the more realistic model of an axisymmetric 
neutron star with a strong dipole magnetic field aligned with the 
rotation axis. In this case acceleration occurs at large distances 



from the surface of the star. The magnitude of the current can be 
determined from this model and is found to be the same as estimated by 
S tarrock (I97I) • In "the case of non- axi symmetric systems the accelera- 
tion is expected to occur nearer the surface of the star. 
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Chapter I 


INTRODUCTION 

1.1 DISCOVERY AND EARLY THEORIES 

The discovery of the pulsars by Professor Hewish and Miss 
Jocelyn Bell in 1967 ranks with the discovery of quasars and 
of the universal microwave background radiation as one of the 
major advances in modern astronomy. 

F.G. Smith. Pu1sars >P. xi 

The discovery of pulsars aroused immediate and intense interest not 
only among astrophysicists and astronomers but the public as well. A 
number of attempts were quickly put forward to explain these remarkable 
objects. Perhaps the most popular with the general public was the so 
called "LGM" (Little Green Men) theory, which suggested that the pulses 
were signals from an advanced extra-terrestrial civilization. This 
theory, however, was quickly discounted; the signals were too regular; 
since the periods were unmodulated they carried no information and it 
was therefore highly unlikely that any little green men were using them 
as communications beacons. It was still possible that pulsars were some 
sort of galactic lighthouses but it was clear that astrophysicists would 
do well to look for a more natural (though less exciting) explanation. 

The most immediately attractive idea was that pulsars were related to 
white dwarf stars and several theories were developed along these lines 
(e.g. Ginzburg et al . 1968; Black, 1969). At the same time, however, 
some work was being done on the possiblility that pulsars were related 
to neutron stars (Gold, 1968; Pacini, 1968). In fact, even before the 
first pulsar was discovered Pacini (1967) had suggested that a 
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magnetized rotating neutron star was responsible for the energy budget 
of the Crab Nebula, The controversy did not last long and was finally 
settled by the discovery of the Vela pulsar (Large, Vaughan and Mills, 
1968) and the Crab pulsar (Staelin and Reifenstein, 1968), which had 
much shorter periods than previously discovered pulsars. The white 
dwarf theories were now running into serious difficulties, which are 
summarized in table 1 below. In the first place, it was clear that 
white dwarf stars could not be rotating with periods any faster than 
approximately 8 sec. (This is the period for which the gravitational 
force equals the centrifugal force at the surface of a white dwarf). 
Vibrational modes of a white dwarf could afso be rejected. Since 
vibrational periods are approximately given by 

P ^ (Gp)-^/2, (1.1.1) 
the expected pulsation period for white dwarfs is of the order of 1 sec, 
which fits reasonably well with the first discovered pulsars but is 
difficult to reconcile with the Crab or Vela periods. In addition, it 
is difficult to understand why only one mode is observed and why the 
mode is so stable*. Furthermore, as a white dwarf ages it cools and 

i 

contracts slightly. Thus the density increases and the vibrational 
period uould be expected to decrease; instead pulsar periods are 
observed to increase, which is what one expects of a rotating system. 
Since white dwarf rotations must be rejected, we are left with neutron 
star rotations. Thus, with the discovery of the Vela and Crab pulsars, 
white dwarf models were no longer tenable and it was clear that neutron 
stars had finally been observed. 
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TABLE 1 


Neutron Star vs. White Dwarf Models 



NS 

WD 

1. Period in range 0.03 to 3.7 s 

V 

X 

2. Period stable to one part in 10^ 

V 

? 

3. Period increases 

V 

X? 

4. No optical photospheric radiation 

V 

X 

5. Two pulsars in supernova remnants 

V 

X 


Having determined that the pulses were due to the rotation of a 
neutron star it was now necessary to develop a more detailed model for 
the emission mechanism. 


In a recent book of the subject of pulsars^ Manchester and Taylor 
(1977) remark: 

One of the least understood aspects of pulsars is the 
mechanism by which rotational energy is converted into pulses 
we observe. Although numerous theoretical models for the 
emission mechanism have been proposed^ no single model has 
been generally accepted. 

‘It is clear that the emission mechanism must be a coherent one. The 
brightness temperature at a given frequency is defined by 

I (v)c^ 

TbCv) = (1.1.2) 

2kv^ 

For typical pulsar parameters this gives brightness temperatures in the 
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range of 10^3 k to as high as Tb == K, For incoherent processes 
this implies particle energies of the order of kTb ^ 10^^ eV, It is 
difficult to imagine an acceleration mechanism which will produce 
particles of such enormous energies and even if such such highly 
energetic particles were produced^ they would radiate most of their 
energy in the frequency band around 5,8 x 10^*^Tb - 5,8 X 10***^ Hz, Such 
energies have never been observed and we may therefore reject incoherent 
processes as the source of the radio emission. There have been many 
suggestions for the coherent mechanism, but none has been completely 
satisfying and I will have little more to say on this • subject in this 
dissertation. 

A second question must also be considered in relation to the radio 
emission. Where is the radiation produced? There are two main schools 
of thought on this question. The first (e,g. Gold, 1969 ) advocates the 
'Might cylinder" model, in which plasma follows field lines out to the 
light cylinder. At the light cylinder, the plasma is highly 
relativistic and radiation is beamed in the forward direction. Light 
cylinder models were refined by F.G, Smith (1971 and 1973 ) but little 
recent theoretical work has been done to develop detailed pulsar models 
with emission at the light cylinder. The alternate model was initially 
presented by Radhakrishnan and Cooke ( 1969 ), In their model, the 
emission region is near the surface of the neutron star, in the region 
above the magnetic polar caps. The radiation is assumed to be beamed 
into a cone (known as the emission cone) and thus acts rather like a 
lighthouse beacon. The relative merits of these two pictures is still a 
subject of some controversy, but a partial summation is given in table 2 



TABLE 2 


Light Cylinder vs. Polar Cap Models 


Light Cy 1 inder Model s Mer i ts 

1. Natural beaming process 

2. Rapid, asymmetric changes in polarization within subpulses 

Light Cvl i nder Model s - Deficiencies 


1. The strength of the magnetic field at the light cylinder 
depends on the pulsar period. Thus pulses from slow 
pulsars might be expected to be very different from 
pulses of fast pulsars. This is not observed. 

2. The emission region is small compared to the light 

cylinder radius. A mechanism must be found continuously 
to supply particles to the emission region while 

maintaining the coherence of the process. 

3. The stability of the pulse shapes indicates that the 

emission takes place in a region of strong magnetic 

fields where the particles co*-rotate w.ith the star. 
This is unlikely near the light cylinder, 

Pol ar Cap Model s Merits 


1. Simple explanation of the stability of even very complex 
pul se prof i 1 es. 

2. The emission region would be expected to be small 
compared to the entire stellar surface thus producing 
pulses with widths of the order of 10® of longitude (as 
observed) . 

3. Strength of the magnetic field in the emission region is 
independent of the period and hence pulse 
characteristics would be expected to be relatively 
insensitive to period. 

Pol ar Cap Model s :i Deficiencies 


1. The simplest polar-cap models predict a pulse width that 
is smaller than observed. 

2. Particles must be accelerated to highly relativistic 
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energies in order to produce and beam the radiation. 
But the charge density in the magnetosphere may be 
expected „to adjust to decrease the acceleration (if 
possible) . 

3. The relation between energy loss and angular momentum 
loss suggests that the primary processes (energy and 
angular momentum loss) affecting the star must occur at 
the light cylinder (C.f. Holloway, 1977). 


The one truly outstanding problem with polai — cap models is the source 
of the acceleration, which is required to produce coherent radio 
emission near the stellar surface. This dissertation is primarily 
devoted to an attempt to deal with that problem. In Chapter II, section 
1 I will present the basic model in more detail, while section 2 will 
deal with an analysis of the problem of particle acceleration in polar'- 
cap models. In Chapter III I present a new approach to the acceleration 
problem using a more realistic magnetic field structure than in previous 
work. Finally, in Chapter IV I discuss the results of this research- 




chapter II 


THE MAGNETOSPHERE PROBLEM 
2.1 THE PCLC AND PCFB MODELS 

Before turning to the main body of this thesis it is necessary to 
define in more detail the salient features of the Stanford pulsar 
models, which form the basis of the current work. The original model 
was developed by P.A. Sturrock (1970, 1 971a, 1971b) and has formed the 
basis of all subsequent development of pulsar models at Stanford. In 
1989, P. Goldreich and W.H. Julian published a paper of fundamental 
importance to the pulsar problem. In this paper they demonstrated that 
"in spite of its intense surface gravity, the star must possess a dense 
magnetosphere." The plasma in the magnetosphere has essentially 
infinite conductivity and hence obeys the "f rosen-in-f 1 ux" condition. 
The magnetic field lines may be viewed as being firmly attached to the 
surface of the neutron star and, as the star rotates, the plasma in the 
magnetosphere is forced to rotate along with it. This cannot, of 
course, be true if the plasma would be forced to move faster than the 
speed of light and hence at the light cylinder the "frozen-in-flux" 
condition requires that magnetic field lines be pulled out and wrapped 
around the star. The basic picture is shown in Figure 1 below. 

The distance, Ry , -is the radius of the "Y-type neutral point" and it 
defines the field line which separates field lines that are closed 
within the co-rotating magnetosphere from lines that are (in some sense) 
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Figure 1: Basic picture of the magnetosphere of an axisymmetric 

rotating magnetic neutron star 


open* and connect to the interstellar medium surrounding the star. In 
the Goldreich-Jul ian model (hereafter referred to as the G-J model) the 
radius is the light cylinder radius defined by 

Ry = = c/n = cP/2tt, (2.1.1) 
where is the angular frequency and P is the pulsar period. It is then 
assumed that particles flow freely along magnetic field lines (i.e. 
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This determines a charge density in the magnetosphere given 


E»B = 0). 
fay 



( 2 . 1 . 2 ) 


The model developed by Sturrock (referred to as the Polar. Cap J-ight 
Cylinder model) ( 1970, 1971a, 1971b) is basically an extension of the work 
of Goldreich and Julian in which the condition that E*B = 0 everywhere 
is relaxed* Specifically, the condition does not apply to the open 
field lines* Thus, on open field lines, particles can be accelerated to 
very large energies* 

The polar cap is defined by the condition that the magnetic field 
line which leaves the edge of the cap be the last closed field line* 
Thus all field lines emanating from the polar^cap region are open field 
lines and particles may be accelerated along these field lines. The 
equation defining a dipole field line is 

sin^e/r = const, (2.1.3) 
The polar cap angle Bp is then defined by 

sin^8p = C2.1.4) 
The rotation of the star induces a potential difference between the 
center of the polar cap and the edge. In the simple case of an aligned 
rotator, the potential on the surface of the star is then given by 


§ = 


- 

2c 


2 ^ 
cos 9 


(2.1.5) 


where B* is the strength of the magnetic field at the pole and R* is the 
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radius of the star. Thenj from equations (2,1,4) and (2,1.5) we find 
that the potential difference from the center of the polar cap to the 
edge is 




2 
2 c 


• 2 ^ 
sxn 0 


fi nl 

2 C 


( 2 , 1 , 6 ) 


For typical pulsar parameters this gives a potential difference of the 
order of 10^^ Volts, Thus in this model we may expect charged particles 
to be rapidly accelerated to highly relativistic energies. In this 
models each polar cap produces two current streams. In particular> if 
> 0 electrons are accelerated from the central portion of the poVar 
cap and ions are accelerated from an annulus around the central area. 
The two zones are referred to as the "electron polar zone" (EPZ) and the 
"ion polar zone" (IPZ) respectively. 


Because the particles follow curved field lines, they emit photons of 
energies 

£(ev) ^ (2,1.7) 

c 

where A is the mass in a.m.u, (A=10"^'^^ for electrons), E is the energy 
of the charged particle (in eV) and Rc is the radius of curvature of the 
field line. As the photons cross magnetic field lines, they "see" a 
changing, transverse magnetic field with which they can interact, 
producing electron-positron pairs (Erber, 1966; Daugherty and Lerche, 
1976), In this model, the pair production process is necessary for the 
production of coherent radio emission. Clearly, for this mechanism to 


work, the 

initial gamma rays 

must 

have 

energ ies 

above the 

pai 1 — 

production 

threshold* 

which in 

turn 

requires that 

the energy 

of the 

initial particles be 

above some threshold 

energy 

(dependent 

on the 
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curvature of the field lines). This then gives a natural explanation of 
the "turn-off" condition for pulsars. This condition has been further 
investigated by Sturrock , Baker and Turk (1976) and has been generalized 
to include radiation reaction and distorted magnetic fields. 

This model (like most polar-cap models) predicted a definite 
relationship between the pulse width and pulsar period given by 

W « (2.1.8) 
As can- be seen, from figure 2, the PCLC model does not fit this 
distribution at all well. In addition, the PCLC model (along with most 
polai — cap models) predicts that the braking index defined by 

n = wcj/(j2 (2.1.9) 
have the value n=3. It is very difficult to determine the braking 
index, but for the Crab pulsar the current best value is n=2.215i.005 
(Groth, 1975). 

This led D.H. Roberts and P.A. Sturrock ( 1972a, 1 972b, 1 973) to modify 
the PCLC model by changing the "Y-type neutral point" from the light 
cylinder radius, R^, to the "force balance radius", which is the 
radius at which the co-rotation velocity is the Keplerian velocity for a 
circular orbit (Roberts and Sturrock, 1972a, 1972b, 1973). 

Rpg = (Gtl) (2.1.10) 
In this model (called the PCFB model), the polar cap angle, Op ,is given 
by 

9p 10’ • 5n.' 1/3 (2.1.11) 
and hence the pulse width is proportional to As can be seen from 
figure 2, the fit is much better. In the region r<Rpg the magnetic 
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field is assumed to vary as r*^ while in the region beyond Rj.g the 
magnetic field varies as r~^. This change in the magnetic field 

structure changes the torque and the braking index is then given by 
n=2.33, which is in better agreement with the observed value for the 
Crab pulsar. 



12 







ORIGINAL PAGE IS 
OF POOR QUALITY 

2.2 SELF-CONSISTENT MAGNETOSPHERES 

The first comment that needs to be made on the subject of self- 
consistent magnetospheres is the limited use* of the term '^self- 
consistenti" To be truly self-consistent, a model of pulsar 
magnetospheres would have to include a) the effects of currents in the 
magnetosphere on the magnetic field structure, b) the effects of 
particle masses on the currents which develop, and c) the effects (via 
the plasma) of radiation produced in the magnetosphere on its structure 
Ce.g., radiation reaction, self absorption, scattering, etc.). That 
detailed a model is well beyond the scope of this thesis. By ''self- 
consistent" we shall mean models which satisfy the appropriate equations 
without inducing large scale changes in the original conditions. 

In 1974 N,J. Holloway published an important paper which illuminated 
some severe problems with the PCLC, PCFB, and similar pulsar models 
(e,g. Hinata, 1973; Hinata and Jackson, 1973). Holloway pointed out 
that there was a fundamental inconsistency in these models. Consider a 
cylindrical "gaussian pill-box" at the polar cap. The flux through the 
bottom surface is zero (space-charge limited flow), the flux through the 
top is nearly zero provided the "pill box" is extended far enough up to 
get it out of the accelerating region, and the flux through the sides is 
given by the co-rotation electric field (vxb). The charge enclosed is 
therefore approximately 

Q = -{n«B/2Tic)TTr^h (2.2.1) 
where rp is the radius of the polar cap and h is the height of the 
"pi 11 -box". However^ the currents from the EPZ and IPZ are expected to 
be comparable and the net charge enclosed should be approximately zero. 
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It is interesting to note, however, that the charge density is 
consistent with the current flow if only one sign of charge is 
accelerated* This suggests that instead of having currents of opposite 
charges both flowing out from the star, we may instead have a current 
loop, with the return current being^ outside of the polar*-cap region* 


It is possible that the large numbers of e'*'“*e' pairs produced may 
adjust their distribution so as to satisfy equation (2.2.1) while the 
two currents flow through the pair plasma. The objection to such a 
model was well stated by Holloway: 


... In the positive particle acceleration zone of such a 
system, there would have to be an electric field which 
accelerated the positive particles to high energies, while 
leaving the negative particles essentially unmoved. a 
situation which, while perhaps not demonstrably impossible, 

(one could postulate a situation in which some form of plasma 
streaming instability counteracted the systematic fields) 
seems at least, implausible. Furthermore, in the regions above 
the accelerating zones, the required coexistence of a 
relativistic, high density, stream of particles, with a static 
corotational charge density of the opposite sign, would seem 
to present great difficulties for this model. 

M.A. Ruderman and P.G*. Sutherland developed ( 1975) a new pulsar model 


which used a very clever idea. A fundamental point of the problem is the 
assumption that the accelerating electric field is zero on the stellar 
surface. Ruderman and Sutherland pointed out that if the work function 
of ions were high enough they could not be removed from the stellar 
surface. Thus, in the case that < 0, so that ions must be removed 
from the central region of the polar cap, a vacuum region will develop 
(called the "polar gap") and a large accelerating electric field will 
form at the surface of the star (and in the entire "gap" region). In 
this model the accelerated particles come from the static breakdoun of 
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the vacuum. Since this model relies on the work function for ions being 
very large, it requires that all pulsars have n*B < 0 (neutron stars 
with n*B > 0 would not accelerate particles and hence would not produce 
coherent radiation). The reason for the large work function for ions is 
that in a strong magnetic field gauss) the ions form long chains 
and the gravitational binding of the chain is large. Recently, Flowers 
and his co-workers (including Sutherland) (Flowers, et al . 1977) have 
recalculated the work function and found that the work function used in 
the Ruderman-Sutherl and pulsar model had been over-estimated by 
approximate! y an order of magnitude. With the new work-function 
estimate, the polar gap does not develop, and the net result is (in 
Sutherland's own wordsM that ^'the model is dead." 

Both F, Curtis Michel and E,A. Jackson have developed pulsar models 
which avoid Holloway's criticisms. They also, unfortunately, produce 
very little acceleration and provide no mechanism for the observed 
radiation. Michel's model (1975) utilizes currents of a single sign 
moving on radial field lines (see Chapter 3, section 1) and thus simply 
matches the G-J charge density. In the simplest form of the model 
(where particle inertia is ignored) no acceleration takes place at all. 
When the particle masses are taken into account there is acceleration 
until the particles become relativistic, at which point the acceleration 
ceases. Even in this case, the acceleration is not sufficient to 
provide a mechanism for the observed radiation. Typical values of the 


^Private communication made to the author at the eighth Texas Symposium 
of Relativistic Astrophysi cs, Boston, MA., Dec. 19.76. 
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relativity parameter y are of the order of 10, compared to 10’° for the 
PCtC model . 

Jackson's model (1978) abandons the requirement of space-charge 
limited flow and substitutes field emission (at T " 0 K) at the surface. 
The current is then related to the accelerating electric field by 

j*= t6.2x10'°E2(ii/?^)’''2/(55+n)]expt-6.8J<107?53/2/E„l (2.2.2) 
where Ej, is in volts/'cm, ^ is the work function in eV, and p. is the Fermi 
energy relative to the bottom of the conduction band. This model also 
features complete current loops, so the requirement of zero net current 
leaving the star can be dropped (since no current at all leaves the 
star). The difficulty is, again, that -there is very little acceleration 
and no reasonable radiation mechanism. 


0? POOR 
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MATHEMATICAL FORMULATION 


3. 1 BASIC EQUATIONS 


3.1.1 The Basic Model and Acceleration Theorem 

The pulsar model presented here is a development of the polar-cap 
models of Sturrock and Roberts and Sturrock (see Chapter II» section 1). 
Fundamental to this model is the fact that currents flou from the polar 
cap along magnetic field lines. That this is true can be demonstrated 
by comparing the gyroradius of the particles uith the radius of the 
polar cap. The gyroradius is given by 

rg = pc/eB £/eB A$/B C3.1,1) 

If we take the maximum that we can get (equation 2*1,6) ue find 




(3. 1.2) 


” 2c 


B 


Near the polar cap B -B^ and the polar cap radius is given by 




'p " J.V2 


(3,1,3) 


1/2 


and hence the ratio of r^/rp is 


r 


r 

P 




^ (3,1,4) 


where Vg 
however, 
magnetic 


is the rotational 
that as r -> 
fiei'd controls the 


velocity of the stellar surface. We note, 
rg Rj^, We can also estimate whether the 
current flow or whether the current flow 
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controls the structure of the magnetic field. For this we simply 
compare the energy density due to the kinetic energy of the particles to 
the energy density of the magnetic field. The particle number density 
is given by 

n = p/e «« HB/(2Tfce) (3.1.5) 
For the energy we again use eA^. Thus the ratio of kinetic to magnetic 


energies is estimated to be 


3 3 

20-^ R* 






(3. 1.6) 


Near the polar cap the ratio is very small and the magnetic 


field controls the particle flou. When r grous to the order of 
however, the ratio approaches unity and in that region we may expect the 


magnetic field to be distorted by the particle flow. 


We assume that the plasma is completely charge separated, which means 
that the pair production process is not taken into account in 
investigating the acceleration mechanism. This treatment would also be 
valid provided the net current due to pairs is small compared to the 
primary current from the polar cap. If the acceleration is large this 
will clearly not be the case and in the region of large pair production 
the model will break down (see Chapter IV). Since the particles are 
tied to magnetic field lines, the current density is proportional to the 
magnetic field strength. Thus, along a field line we may write 

j(s) = j(0)[B(s)/B(0)I (3.1,7) 
where s Is a co-ordinate along the field line and s=0 refers to the 
surface of the neutron star. Since the particles are relativistic (as 
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Mill be demonstrated in section 2), the charge density is given by 

U B 

p (3.1.8) 

c 

The G-J charge density (equation 2.1.2), uhich is required for ^*^=0 
(i.e. no acceleration), is proportional to fl*B. Thus, if the angle 
betueen jl and ^ is constant along a field line then equation (3.1.8) is 
compatible with (2.1.2) (in the non-rel ativistic limit) and it is 
therefore possible to have steady current flow with no acceleration. 
If, however, field lines curve then acceleration (or deceleration) must 
take place. This theorem has been independently derived by Arons, 
Fawley and Scharlemann (1978) by transforming to a rotating reference 
frame^. In the frame rotating with the star, there is an electric field 
perpendicular to the magnetic field, given by (in the non-relativistic 
limit) 

Erot - KSl X r) X Bl/c (3.1.9) 

If we assume that B is approximately curl free (i.e* the magnetic field 
of the star is much larger than the field generated by currents in the 
magnetosphere), then the charge density of ^rot is given by the 6-J 
charge density (equation 2.1.2). Ue may then divide the electric field 
into two parts, the rotational part given fay equation (3.1.9) and the 
non-rotational part uhich may accelerate part.icles. Thus only the 
difference between the rotational charge density and the true charge 


^Tademaru (197<1) proved a restricted version of the theorem, too. He 
showed that for an axisymmetric rotating system with a polai — cap region 
bounded by radial (i.e. monopole like) field lines, the component of E 
parallel to the magnetic field lines must be exactly zero. 
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density is a source of acceleration (Michele 1975). When the angle 
between SI and B is constant we can find a current flow such that the 
noM'-rotational charge density is zero everywhere. When the angle is not 
constants however, that is no longer possible and hence a non'-rotational 
electric field must develop. We are thus motivated to look at a pulsar 
model in which the field lines are curved. The simplest physically 
realistic example is a pure dipole field. 


3.1.2 Dipole Co-ordinates 

In order to study the dipole field case we first introduce a co*~ 
ordinate system based on the dipole field lines. The potential of a 
magnetic dipole oriented along the 2 axis is given by $ 0 : cos 8 /r^. 
The equation of a field line is sin^6/r = const- Thus for our co- 
ordinates we may take 

f = Vcos0/r and 7) = sin8/\Tr (3.1.10) 
The third co-ordinate is the azimuthal angle but we shall usually 
assume azimuthal symmetry and thus reduce the problem to two dimensions. 
The Laplacian for dipole co-ordinates can be written as 




1 

T1 



C3. 1. 11) 


Unfortunately, the values of r and 0 cannot be expressed explicitly nn 
terms of ^ and 7), so the Laplacian cannot be written simply in terms of 
the dipole co-ordinates. 


The polar-cap region is defined as the region of open field lines. 
The bounding field line is determined by the value of Ry» the "Y-type 
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neutral point" (see equation 2.1.11). Thus we can define a maximum 


value for the co-ordinate 7 } by 




(3.1.12) 


Since we are assuming a charge separated plasma, the current density is 


j (r) = p(r) V (r) 


(3.1.13) 


Since the particles are tied to field lines, j can be separated into two 
components^ the motion along the field line and the motion due to the 
rotation of the star. However^ motion due to rotation is small compared 


to that along field lines if r << (the light cylinder)* t4e 


shal 1 


therefore neglect the rotation part of the current d^nsity^ in which 
case equation (3.1,13) becomes a scalar equation. Combining it with 
equation (3.1.7) we then have 




j»(T|) lg(g.'ri)| 

vCl.Tl) iB(|*,T|)| 


(3.1.14) 


where is the value of ( on the surface of the star (f^ ^ snd v 
is the velocity of the particles. Since the field is dipolar# we can 
write as 


B(E,7l) = ^ (3- - f ^ 

r ' ' 


(3.1.15) 


where p, is the dipole moment and r is implicitly determined by $ and 7 ). 
Combining equati ons_(3. 1 . 14) and (3.1.15) we then have 


, ^ j,(u) /»*\ri 

P(I.H) * y(§,1l) ^ 


(3.1.16) 
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For many purposes we can ignore the dependence entirely. Since 

r7)^ i 1 we can make the expansion 




3 2\ 



and even at the extreme value oi rD^ = 1 C6 = tt/2) 
no more than 25%. 


(3.1.17) 
ue make an error of 


3.1.3 Boundary Conditions 

Ue must now consider the boundary conditions appropriate to the 
problem. One boundary condition is clear: the bounding closed field 
line satisfies the condition ~ ^ hence the potential along the 
field line must be a constant, which we may take to be zero. In terms 
of dipole co-ordinates we therefore have 

= 0. (3.1.18) 
(Unless otherwise stated, # refers to the electric potential.) The 
surface of a neutron star is a good conductor and therefore the Lorentz 
force on a charged particle on the rotating surface must be zero. Hence 
the electric field parallel to the stellar surface must be the rotation 
electric field (equation 3.1.9). If we assume that ^ is the gradient of 
a potential with azimuthal symmetry, we can write 


“ Pj^(cos 9) 


(3.1.19) 


for the magnetic potential, where li is the moment associated with the 
/-pole. Thus $, the electric potential can be determined on the 
conductor surface by 

i = -R*XUnxr)x17$jj]de (3.1.20) 
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If the rotation axis is aligned with the symmetry axis of the magnetic 
field this is easy to integrate. Since the velocity of the surface is 
entirely in the ^ direction we need only consider the radial derivative 
of the magnetic potential. We then find that the electric potential is 
given by 




E 


2^ \^j — 

-t=l ^ 


(3.1.21) 


In the case of a pure dipole magnetic field this becomes 

§ = (cos^9 cos^9p)/2c. (3.1.22) 

In terms of the dipole co-ordinates we may rewrite this as 


$ 


E 


n 


r3 

* %ax 
2c 



(3.1.23) 


The remaining boundary condition is more difficult to determine. 

Typically, the remaining boundary condition would be at infinity (r=<»). 

In the case of dipole field lines, however, the field lines do not 

extend to infinity. Nevertheless, at r=«> it is true that ^=0 but | is 

also zero at 0=ti/2. The dipole field structure must breakdown at the 

light cylinder (or perhaps at some Ry < Rj^) and the "frozen-in-flux" 

condition also breaks down when E > B, a condition we may again expect 

near the "Y-type neutral point." This indicates that we should be 

looking for a boundary condition which applies at the point r=Ry. If 

the "frozen-in-flux" condition breaks down at Ry it is no longer 

reasonable to assume Ei - -v/c x B. In addition, since B " r'® up to 

^ ^ 

the radius Ry, B(RyX<B*. Therefore (as pointed out at the beginning of 
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this chapter), the gyroradius 1s comparable to Ry and particles can 
cross field lines, effectively shorting out the c1rcu.it. Hence we may 
consider the approximation that " 0 at Ry and therefore the electric 
potential is approximately constant at Ry. Since this surface 
intersects with the bounding field lines, on which the potential is zero 
we must have 

§(Ry,0) « 0 (3.1.24) 
which then provides a complete set of Oirichlet boundary conditions. 

There is one additional condition, however, that we may be able to 
apply. At the surface of the star (where particle flow is not 
relativistic) a region of charge may develop which then decreases the 
electric field normal to the surface. This phenomenon is well known 
from vacuum tube technology and is known as "space-charge limited flow". 
If the accelerating field were non-zero at the surface we would expect 
the current flow to be increased (a stronger electric field would pull 
out more particles) and hence the charge density would increase. Steady 
flow is achieved when the current flow is just sufficient to keep the 
accelerating field zero at the point of particle emission. Space-charge 
limited flow is an important feature of the PCLC and PCPB models and has 
also been invoked by many other investigators (though not all, e.g. 
Jackson, 1976) . Thus as an additional constraint on the problem we will 
consider the case of space-charge limited flow, which implies that the 
accelerating electric field at the surface must be zero. 

U = 0 (3.1.25) 
We h0J4 have an overdetermined problem and we are no longer free to 
choose the Initial current density at the surface of the star. The 
equations and boundary conditions are collected in table 3. 
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TABLE 3 


Equations and Boundary Conditions 








! 3 2 

p - f w 


( 3 . 1 . 16 ) 


1 - ^rrt® 



iliLii. a 


1) 9T) I ' 




w) = 0 


§(Ry,0) = 


,, ^ ^ %ax / _ a . 2 



3.2 ANALYTIC AND NUMERICAL SOLUTIONS 
3.2.1 Th^ Nj^p- linear Prob1 em 

We uil? first show that the partfcies from the polar cap are quickly 
accelerated to relativistic velocities. In this analysis we will 
consider a one-^dimensional model (that assume the divergence of E 
across field lines can be neglected with respect to the divergence along 
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the field lines) and a small angle approximation. Poisson's equation is 
non-linear in this case because the velocity in equation (3,1,16) is 
dependent on the potential. Rather than work with the potential we 
shall work directly with the particle energy. 

e = rmoc^ = moc^ - e$ + const. (3.2.1) 


so 


= -(mocVe)V^r 

The velocity can be expressed in terms of y as 


V = 


/ 2 xV2 

(v - 1) 


(3.2.2) 


(3.2.3) 


and hence Poisson's equation reduces to 




V 


^ - 1 ^ 
1 - i 


-, 1/2 




(3.2.4) 


We now assume that (at least until the motion becomes relativistic) the 
divergence of the electric field along the magnetic field lines is much 
greater than the divergence of the electric field perpendicular to the 
magnetic field lines. Thus t) derivatives can be ignored. We also make 
a small angle approximation and ignore terms. Equation (3.2.4) then 
becomes 



(3.2.5) 


where the dipole nature of the field lines is reflected by the r^ 
dependence. Equation (3.2.5) cannot be solved analytically and so 
additional approximations must be made. The necessary approximation is 
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to assume a cylindrical model for the structure of the magnetic field 
lines near the surface. Thus equation (3.2^4) becomes 


dr2 -VJ (y"-!)'-/" 


(3.2.6) 


This can be integrated once by multiplying both sides of Poisson's 
equation by dy/dr giving us 


dY 

dr 





(3.2.7) 


This may then be integrated in terms of elliptic integrals. We first 
make the substitution 

cosh y = 7 . (3.2.8) 

Then the solution to equation (3.2.7) is given (implicitly) by 


, \l/2 

2(sinh y) cosh y 

1 + sinh y 







r " 

jl/2 



8jtej^ 


1 

1 

1 

3 

0 
o 

1 

_n/% 



(3.2.9) 


where F and E are elliptic integrals of the first and second kind 
respectively, and 


^ /I - sinh y\ 

a = arccos I ^ ; — — 1 = arccos 
\1 + sinh yj 


(3.2.10) 


1 + 


7y^-i 


We must remember, however, that this equation is valid only for r - R* 
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so that the divergence of the field lines could be ignored* Rather than 
working with the full analytic result/ it is easier to consider the 
additional approximation that y ^ 1 (i.e* non-rel ati vistic motion)* We 
therefore write 

y = 1 + 8 (3,2.11) 

and keep only terms to first order in S. Equation (3.2.7) then becomes 


Si 

dr 





(3-2.12) 


Integrating, we find 


i a3A 
3 ® 


,5A 




11/2 


mQC' 


R 


3/2 


* 



1 


(3.2.13) 


and since we know r we can further approximate the right hand side 

to get 


3 


A/if 


8jrej, 


P -/2 


V j 


(3.2. 14) 


where r = R^+ Sr and Sr << R*. Finally, we can solve for 8 to get 


3^/3 

A/3 


Jtej 


- 2/3 


* 


V 


A/3 


(3.2. 15) 


Typical estimates for the current are of the order of 10'^ esu/cm^-sec. 
Thus, we write equation (3.2.15) in the form 

S .56j2/3(Sr)^^.3 (3.2.16) 

where we have assumed the particles are electrons and where 
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ji 2 - jit ^ 10*'^^ and Sr is measured in centimeters. In the non- 
relativlstic limit, S - and hence when S - i the particles have 
velocity \ - c and the non-rel ativistic approximation is no longer 
valid. It is clear from equation (3.2.16) that the particle becomes 
relativistic within a few centimeters to a few meters (depending on the 
magitude of j*). We can also estimate the distance over which the 
divergence of the field lines becomes significant. From the approximate 
solution (equation 3.2.15) we find that 

— (3.2.17) 
dr 3 ^ ^ 

For typical pulsar parameters, Poisson's equation now can be written 
approximately as 

0.3 - 2 X lO"^ (3.2.18) 
dr 

where the second term on the right hand side had previously been 
ignored. The second term becomes comparable to the first when 
8r “ 7X10“^ cm. Thus it is valid to ignore the divergence of the field 
lines in the non-rel at i vistic limit. 

Having confirmed that the particles quickly become relativistic, it 
is now possible to deal with a linear partial differential equation by 
simply replacing v with c in Poisson's equation, which now becomes (in 
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3.2,2 One Dimensional Solution 

As a first approach to the linear Poisson equation we again consider 
a one-’dimensional approximation; that is we wish to consjder the case 
where the divergence of JE perpendicular to the field lines is small 
compared to the divergence of ^ along field lines. 

jv-Exl << !v-E„| (3.2.20) 
We recall (equation 2.1.2) that the charge density for the case of ^ii= 0 
(no acceleration) is proportional to This suggests that we look 
at the perpendicular rotator^ for which case - 0. In this case, a 
significant departure from = 0 must indicate an accelerating electric 
field. We note that the one-dimensional model must be treated as an 
initial value problem rather than a boundary value problem. Before 
proceeding we also note that because the co-rotation charge density is 
approximately zero near the polar cap, any charge density due to the 
emitted current supports the accelerating field rather than the co- 
rotation field. This would tend to indicate that the acceleration 
probably occurs closer to the star surface in this case than is the 
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general rule. Equation (3.2.19) now reduces- to 








1 - I 


n 1/2 


(3.2.21) 


Simplifying this equation we then find 





r: 


■ I" \ 


1/2 


r2(cos (1 


3 2 

IT ^ '’'1 ) (3,2.22) 


This can be integrated once exactly by integrating along a field line. 
Since 'b^y'd^ = 0 on the surface, we have 


1 

I 9| 


rI 


1/2 




j r^(cos 0)^^^ (1 - I r 7 ]^) d? ’ (3.2.23) 


I* 


From the definitions of f and 7), we may write the differentials 


d| = - 


Vcos 9 ^ sin 0 

V (Jj. 


d© 


2rv/cos 0 


(3.2.24) 


dT| = - 


sin 9 , . cos 0 

dr + — d© 


2r 


3/2 


1/2 


We know di5 = 0 along a field line, and we can therefore express d9 in 
terms of dr to give us 

L . . 


d? = - 


r^(.cos 9)3/2 


dr 


(3.2.25) 
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This is not, - however, a particularly useful form. Since we cannot 
expect the one-dimensional approximation to be valid for large r (and 
hence large 6) it is more useful to return to equation (3.2.28) and 
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expand to second order in 6 (i.e. r7}^). We then have 


1 ^ 

I 




(3.2.30) 


It is now easy to do the next integration to order 0^. We note that 

de = -dr/r2 + (O'*) (3.2.31) 

and 



Thus 



To zero order in we then have 


(3.2.32) 


(3.2.33) 


E 


+ 


dr 



(3.2.34) 


It is interesting to note that the electric field has a maximum at 
r ~ 1.5R*, which suggests that the maximum particle acceleration may 
occur at distances of the order of a stellar radius above the polar cap 
rather than at distances of the order of the polar cap radius as in most 
previous models. We further note that this form for is due to the 
dipole nature of the field lines and is not found in cases where the 
field lines are not curved. 
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3.2,3 Numerical Sol uti on to the non-1 inear 1-D prob1 em 

We now turn to a numerical solution to the non-linear one-dimensional 
model. The equation we wish to solve numerically is 







/ R. 


.3 r 


'* 


T 3 ^2- 

1 - Ij: r 71 


, 1/2 


moC- 


3 2 

.1 - F M J 


(Y^ - 1) 


1/2 


(3.2.35) 


This equation can be solved easily using standard differential equation 
solving programs Cin this case the program ODE^ developed by Shampine 
and Gordon - see Appendix B - was used) provided accurate starting 
values can be determined. The infinity at must be avoided by 
starting the integration at a position slightly above the stellar 
surface, where y~1+8. The results of the first section are used to 
determine the starting values. The second order equation must first be 
decomposed into a pair of first order coupled equations. 
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We now rewrite the system in numerical terms. 


i 

' df“ ^ ^2 








„3 .2 


5 Y. 




21^/2 

- -T5 r T1 J 


Y (3.2.37) 
_2 

r 


where jri = j*xl 0~’2 and R^g= R^xlO'®. To simulate the fact that no 
acceleration takes place on the bounding field line ('>)=7)maxJ we shall 
replace jiz by so that only a portion of the current 
causes acceleration. 

The results of the integration are shown in figures 3 - 6 We note that 
the exact numerical results are completely consistent with the analytic 
approximate results. The initial behavior of 8 with respect to Sr is 
correct and there is> indeed^ a maximum in the electric field at 
approximately 1 . 5R*. 
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It is certainly 


3-2*4 Tmo Dimensional So1 ution 

Me must now consider the two dimensional problem* 
true that the time and dependence can be eliminated in the axi- 
symmetric case, but it is probably a good approximation to use a two 
dimensional model even in the non^axisymmetric case. The two- 

dimensional problem is, however, much more difficult to work with than 
the one-dimensional case. In this section we will consider two 
different approaches to the problem. In the first method we shall 
assume a form for the transverse behavior of the electric potential and 
use a perturbation expansion in the co-ordinate t). The alternate 
approach is to assume that the current has the same 7 ) dependence as the 
potential and use a separation of variables technique* 

3. 2. 4*1 Perturbation Method 

Me first analyse the two dimensional problem from the perturbation 
expansion approach. We know that the potential for the aligned rotator 
must be an even function of 7) so we expand it in a power series in 

= Z ai(f)7)2i/7)2i (3.2*38) 

max 

We also knou that at the surface of the star the potential has the form 
) which thus forms the boundary condition that the 
coeff iecients of all terms of order 7 }'* and higher must go to zero at the 
surface of the star, while ai must go to 1 (we may set ao=1 without loss 
of generality). Finally, the condition that ^ = 0 at TJmax requires 
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00 

^ aid) = 1 (3.2.39) 

i=l 


The basic assumption ue need for a perturbation expansion is that 
[a,* I < I a,*-i I . Even with this assumption the problem involves a great 
deal of messy algebra. In working with the perturbation expansion, we 
nust also expand 3 ^ as a power series in 7 )^. 

Ue now look at the problem to zeroth order in 7 ). The potential must 
be expressed to first order and the current to zeroth order in 7 j^. Ue 
therefore set 

i = $(?)n-7)^/'7)2 ) (3.2.40) 

max 

Substituting equation (3*2.40) into equation (3.2.19) we find (to order 

I 

7?^) 









(3*2.41) 


The zeroth order equation is then 


1 a /l 

K U 


“5" 

r I 


3.„2 

^ ^Max 




0 


(^y 


(3.2.42) 


A particular solution to the inhomogeneous equation (3.2.42) is 




0 „3 ...2 




'Max 


(3.2.43) 


42 “ 



The homogeneous equation may be written as 



To solve this equation we make the substitutions 

u = arid = u^yCu) 

Equation (3.2.44) then becomes 


(3.2.44) 


(3.2.45) 


d^Y ^ 1 ^ 

^ 2 u du 
du 



(3.2.46) 


This is the modified Bessels's equation of order 2 and parameter 4/7}n,ax- 
The solution to equation (3.2.42) is then given by 






/it^y 

l^Max / 









2 

Max 


(3.2.47) 


The boundary conditions at r = Ry and r = R* now determine the values of 
the constants. If Ry >> R,^ it is simpler to make the approximation that 
I - 0 at the outer boundary. Then the I 2 term goes to zero and the 
remaining two terms must cancel. As x 0, 2/x^ and we 
therefore find that 


^2 ^ 

8 ”*^Max ^ c 


r3 

* ''Max 


0 


(3.2.48) 


and hence 


8rtJo rJ 

c 


(3.2.49) 
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At ? = 1/R* we require 


ilfC?*) = 


2 

%ax 


2c 


(3.2.50) 


Since ’■'^/T)max))L we can use the asymptotic expansions of I 2 and Kg 
to simplify equation (3.2.47) at the stellar surface. 


X 

2 


(2<) 


T72 


and 




[kr 


e"^ (3.2.51) 


where 


? = 4^<^2/7J„ax 


(3.2.52) 


Combining equations (3.2.49), (3.2.51), and (3.2.51) we then find 
■Q ^ ^ ^ \l/2 

c.. 


n p p / Tf _r 

=1 = - [— V* ”V*Vx - ®”^0 ( 2 c?) ^ . 


(3.2.53) 


[: 


X I — e“^* 


Finally, we require space-charge limited flow, which determines jo. We 
write d$/df in terms of J (as defined by equation 3.2.52) 


I = ^ IC "Max '^ax <3.2. 

^''Max ^'Wx 


54) 


We also note that 


^ [c^l2(o] = C^Ii(C) and' ^ |c%(C)] = - cV« 


55) 


Hence, at we have 



ind, combining equations (3*2*56) and (3,2,49) with the asymptotic 
expansions (equation 3.2^51)p ue get 


= « e 


-2C 


8 «‘ 


2 . 




— R* e"^^* 


(3.2.57) 


’inally, combining (3.2.57) and (3.2.53), we have 

1/2 


•^0 


(gir) » 4a. 

16:t - R, 


Q. B 


* 


2it 


(3.2.58) 


which is the value of j* originally estimated by Sturrock (1971) and the 
value of j which gives the G-J charge density at the surface. 


One of the original motives for analyzing dipole field lines was to 
get large acceleration near the stellar surface. However, the zeroth 
order solution for the aligned rotator does not produce significant 
acceleration. Actually, this is to be expected, since at this order we 
have only included corrections of order 9^ and ue have simply recovered 
the G-J charged density in a region where the field lines are very 
nearly straight (see the theorem described in section 1.1 of this 
chapter). Ue must therefore extend the analysis to second order in 7) in 
order to determine what acceleration (if any) is produced near the 
surface. Ue set 

I . Kl) [l - * 02(5) - Cgd) J 

(3.2.59) 

j* = Jo 

and we note that az is a function of f but b is a constant. 
Substituting into equation (3.2.19) we now find the zeroth order 
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(l Sf 






2 2 


2 

I 11 




'Max-^ 


C3.2.60) 


while the second order equation is 


!*1 


(«) 




'Max 


^2 iPilf ^2 li!f . 2 li|r ^2 

ri^2 



1 

^^2 

(l 

-.8 \| 


5^ *5 C J 

2 

%ax 

^2 2 
^ ’'^Max 


''Max / / 

(3.2.61) 


kTti r: 


’* 


3 1 3 


b 

2 

"^Max- 


We note that the zeroth order equation has been modified from equation 
C3.2.42) but we are assumming that a^ 1 so to solve these equations 
we first set a£ = 0 in equation (3.2.60) (thus recovering the original 
zeroth order solutions equation (3.2.47)) and then use that solution in 
equation (3.2.61) to eliminate the $ dependence. To solve equation 
(3.2.61), we first define a new variable 

B = (3.2.62) 
With this definition and using equation (3.2.60), equation (3.2.61) can 
be rewritten as 



A_-_2L_/3 
2 2 I? 

''^Max ''^Max ^ 




4 

c 



3 

F 




b 

2 

”^Max 


(3.2.63) 
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4e again change the independent variable to and find 


o du ■ ^ ■ 


du 


"Max 




3 I 3 p b + 1 

^ ^Max 


Tl^ + 
'Max _2 




‘Max 


(3*2,64) 


The boundary condition for B at the surface is clearly B(u^)=0* (Since 
ae=0 at the surface.) At the outer surface the condition becomes 
Bu^ -> 0 as u 0* This is true as long as az does not blou up at 
While this leave's B undetermined on the outer boundary, we shall require 
the stronger condition B=0 at u==0. The condition of space^charge 
limited flow requires 




(pi) 


- 6 + f ^ 


_ R _ 

“ P 2 du 


= 0 


(3.2.65) 


1* 


u„ 


Since B itself is also 0 at u = u* we then have dB/du = 0 at u*. 


We can immediately determine the constant "b" in equation (3.2.59).. 
At u* equation (3.2.64) reduces to 


(b + l) + ' ' ' = 0 


T|. 


'Max 


and thus "b" is given by 


b = 


H-(|J 


‘‘XJq 


- 1 


(3.2.66) 


(3.2.67) 


^Max 

Substituting from equations (3.2.58) and (3.2.50), we immediately find 
that b=0. Thus there is no second order correction to the current. To 
solve equation (3.2.64) we use a Green function approach. The Green 
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'unction for this equation (Bessel's equation) is given by (see Appendix 
, for details) 


. IT 

G(U/U') = — Ja(Xu^) Ja(Au-) + aNa(Au>) 
2a •- • 


(3.2.68) 


there 


md 


a = -/28 


(3.2.69) 


X ~ S/7)max (3.2.70) 

a = -Ja(XUj^)/Na(Xu*) (3.2.71) 

Thus the solution to equation (3.2.64)> with B=0 on the boundaries, is 
given by 


u 


P(u) = - 


“’'Jo 4.x 


* 


-j G(u,u 

0 


0 


3 , 3 „ ,1 , 

^ ^ R* + ^ ^ 

^ '*^Max "^Max 


du' (3.2.72) 


Explicitly substituting (3.2.68) into this equation we find the formal 
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solution 
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3. 2. 4. 2 Separation of Variables Method 

In this method ue look for solutions to equation (3.2.19) in which 
fehe 7) dependence of § is the same as that of j^. Ue further assume that 
I and can be written as 

§(|,7)) = 'fU)H(T)) (3.2.74) 


j*(7)) = joH(7)) 


(3.2.75) 


4e note, however, that equation (3.2.19) is clearly not separable in 
these co-ordinates and ue must therefore approximate it in order to 
render it separable. We write equation (3.2.19) to lowest order in 7 ) to 
get 


^ ll. i. I- /-n 
^ / 11 ^ V ^ j 




■ * 


(3.2.76) 


We now substitute equations (3.2.74) and (3.2.75) into (3.2.76). 




« / // 1 / ‘'O 3 


c * 


(3.2.77) 


Dividing by and regrouping terms gives us 


I. 


^0 * - ./ILh- IL\ 


(3.2.78) 


The left hand side of equation (3.2.78) is a function of ^ only, while 
the right hand side is a function of 7} only. We therefore have the two 
equations 


,/ 2. ^^^0 

|\lf -lit -atlr= 


(3.2.79) 
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and 


+ — h' + = 0 (3.2.80) 

T) 

The solution to equation (3.2.80) is then given by the Jo Bessel 
function. 

H(t)) = Jo(«7)) (3.2.81) 

Equation (3.2.79) will be recognized as identical in form to equation 
(3.2.42) and ue can therefore immediately write down the solution to 
equation (3.2.79) as 

Hi) = ^ ^ (aVr)+ ^4 (3.2.82) 

^ ' ' ' a c 

The boundary condition on closed field lines requires 

Jo(®7)max^ ^ (3.2.83) 

Thus the separation constant a is determined by the zeros of Jo. 


The only remaining task is to evaluate constants. The potential may 
now be expressed in the following form: 


i=l 


■=11 5^2 (’‘l * ‘=2x5'^ i\ ^ 

\ 'Max / y 'Max y 


‘‘”^1 3 2 


(3.2.84) 


^0 (Xi Vv,;) 


't/ll 

where Xi is the i zero of Jo* Following the same lines as in the 
perturbation method (see equations (3.2.49) - (3.2.58)) and using the 
orthogonality properties of the Bessel functions, we find that the 
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32n 

2 ^ 
32 Jt 

2 

Xi c 


"'^Max 


The requirement of space-charge limited flow then determines the ji's. 




and the total current is given by 


■ ‘ 
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chapter IV 


DISCUSSION AND CONCLUSION 

1.1 THE ACCELERATION PROBLEM 
k 1 . 1 The charge separated model 

As is clear from the preceding chapter^ the exact mechanism by which 
charged particles are accelerated is still a problem. Fundamental 1 y 
this is due to our lack of knowledge about conditions in the outer 
magnetosphere (near the- light cylinder, or perhaps near * We know 
that at the light cylinder, the field lines must ^'slip" through the 
plasma or be highly distorted* If this were not true no EflF could be 
generated. An additional complication comes from the requirement that 
acceleration take place near the surface of the star (this is required 
in polar-'cap emission models) and that of space-charge limited flow. 
The space-'charge limited flow produces zero acceleration at the surface 
and the slow change of the angle between H and ^ means that only a very 
small accelerating electric field can develop near the surface. In the 
context of a model with dipole field lines, it seems that we must either 
abandon space-charge limited flow or the requirement that acceleration 
take place near the surface. Larger acceleration will occur near the 
stellar surface, however, if the magnetic field lines near the surface 
curve more rapidly than in the dipole case. 

The latter of these two choices is the more attractive. Referring 
again to figure 2 in Chapter II we note that the data indicate that the 
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polar cap is larger than estimated in the PCLC models* In the PCFB 
Tiodel, the polar cap was made large by moving the point Ry from the 
light cylinder to the force balance radius* An alternate approach, 
however, would be to move the radiation producing region out from the 
stellar surface* The fanning out of the magnetic field lines would then 
produce a wider beam. If this is the correct solution, the data suggest 
that the emission region is at a significantly larger distance from the 
surface than previously suggested in polar-cap models* 

A'simple linear least-squares fit to the data shown in figure 2 gives 
the following formula for the pulse width for a given period: 

logCW) == 1*43 + 0*70xlog(P) (4.1*1) 
where W is measured in milliseconds and P in seconds. It is impossible 
to confirm the slope, however, since there is a range of slopes which 
give a good fit provided the constant term is chosen right* In 
comparing the fit for a slope of i (PCLC) with the fit for a slope of 
2/3 (PCFB), we find there is no significant improvement. The least- 
squares fit for a slope of i gives a constant term of 1*40* The width 
of the emission cone (in seconds) for emission occurring at a distance 
re from the center is given by 



The 1 east-squares fit then suggests that the emission region is located 
approximately at r© ^ 50R^. 

Another difficulty in analyzing the acceleration mechanism is due to 
the nature of the curved field lines. We know that we must have curved 
field lines in order to get acceleration, but the most realistic case 
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(dipole field) is very difficult to work with because Laplace's equation 
is not separable* Small angle approximations are possible (as was done 
in the final section of Chapter III)> but the method is not valid for 
large distances, which appears to be an important region for the physics 
of the problem* A possible solution to this difficulty is to use some 
co-ordinate system which is separable and still mimics the general 
morphology of the dipole field lines* One such co-ordinate system is 
that of ^'toroidar' co-ordinates. Laplace's equation is separable 
(although the Helmholz equation is not) and the 'Afield lines" in this 
system would curve in the same direction as the dipole field lines, 
meeting the equatorial plane at right angles to the plane (as do dipole 
field lines). The radius of curvature is, however, much shorter than 
the dipole case which would probably produce excessively high 
acceleration. The method might be useful, nevertheless, as an 
indication of when acceleration can be expected and perhaps might serve 
as an upper bound on how much acceleration can be expected. 

4.1*2 The effects of pair production 

As noted at the beginning of Chapter III, it is important that the 
plasma be charge separated. Thus, either pair production must not take 
place, or the effects of pair production must be negligible in order for 
the equations to be valid* If pair production does take place, there 
are several effects on the equations. First, the current is no longer 
tied to the magnetic field lines, since a particle may produce a gamma 
ray via the curvature radiation mechanism, which then crosses field 
lines until it produces an electron-positron pair on a new field line. 


55 - 



ORIGINAL PAGE IS 
OF POOR QUALrry 


Thus current can be transported across field lines even if individual 
articles are firmly attached to the field lines. A second effect is 
hat the charge density may no longer be simply related to the current 
ensity. Indeed, once pair production takes place, the number of pairs 
xpected is large, so in order for the charge-separated equations to 
till be valid, the velocity difference betueen electrons and positrons 
lUst be extremely small. Since, however, the electrons and positrons 
lill be accelerated in opposite directions (if the accelerating field 
is non-zero), it is unreasonable to expect the velocity difference to 
remain small even if the initial difference is negligible. We would 
then expect the electron-positron plasma to distribute itself so as to 
shield the particles from the accelerating field. Thus, once pairs are 
produced, the acceleration will be sharply reduced. This might not 
affect the radio emission process (in fact it. may be the source of the 
radio emission mechanism, e.g. the beam-plasma instability postulated in 
the Ruderman-Sutherl and model), but the decreased acceleration makes it 
difficult to account for the high-energy gamma rays that have been 
observed from the Crab pulsar (Ogelman, et al . 1976). 

4.1.3 The effects of particl e inertia 

One of the results of the acceleration theorem presented in Chapter 
III is to suggest that the effect of particle mass is even more 
important that previously realized. It has long been suggested that the 
accelerated particles would tend to wrap up the magnetic field lines 
when inertial effects became important compared to electro-magnetic 
effects. Since it is now clear that such distortions of the magnetic 
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field are accompanied by accelerating electric fields, it is important 
to consider the particle masses in explaining the acceleration 
mechanism* The analysis of the problem with particle masses included is 
extremely complicated and will probably require detailed computer 
modeling in order to solve the equations, 

^•1*4 Return currents ; Is a pulsar charged ? 

In the polai — cap region, particles of only one sign are accelerated, 
and the same particles are accelerated from both poles. It is clear 
that the accelerated particles cannot be accelerated indefinitely and 
removed from the pulsar, or the pulsar would quickly develop a charge 
sufficient to turn off the current. Either particles of both signs must 
be removed from the surface (an unlikely situation except in the 
orthogonal rotator case, see Chapter II - Holloway's analysis) or there 
must be a deceleration region in which the particles are slowed down and 
transported to other field lines in the closed field region and returned 
to the star. One would suspect that particles cross field lines when 
1 e1 > Ib], If the star has a net charge, Q, and we ignore the co- 
rotation E field, the particles would be expected to cross field lines 
when 



If we then demand that particles cross field lines at we can solve 
for H to find 

Q ~ nB*R^/c (4. 1.4) 

In fact, this large a charge on the star itself would be sufficient to 
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stop the acceleration of particles entirely and ultimately leave a 
/acuum around the pulsar (Jackson, 1976). Presumably then a smaller net 
charge appropriately distributed in the magnetosphere, along with the 
rotational electric field is sufficient to make the accelerated 
particles cross field lines. By dealing uith a vacuum "magnetosphere", 
Jackson (1978) has estimated the charge on the star to be a third of the 
above value. In fact, if the charge imbalance is distributed near the 
"Y-type neutral point" rather than on the star itself (as is likely), 
the charge imbalance may be significantly less than either estimate. In 
any case, it is quite likely that a charge imbalance is present and that 
it'forces the accelerated particles to cross magnetic field lines and 
return to the neutron star along the closed field' lines. 

4.2 IHE RADIATION PROBLEM 

While the radio emission of pulsars seems to fit a fairly simple 
conceptual model (although the details of the emission mechanism are not 
well understood), the other forms of radiation present a bewildering 
variety of properties. Two pulsars, the Crab and the Vela pulsars, have 
been detected optically. The Crab pulsar has also been detected in the 
infrared and the x-ray region (Cocke, et al . 1969; Wallace, et al . 
1977). Gamma radiation has been detected from four pulsars 
(PSR-0532[the CrabI, PSR-0833[the Vela], PSR-1747, and 
PSR-1818) (Thompson, et al . 1975; Buccheri, et al . 1976; Ogelman, et al . 
1976). In the case of the Crab pulsar the gamma ray energies are known 
to exceed 1 Gev. All the pulsed radiation of the Crab pulsar (except to 
extremely high-energy gamma rays " 10’^ eV) is in phase with the 


58 - 



>rincipal radio pulses. 


In the case of the Vela pu1sar> there is only 


)ne radio pu1se> but tuo optical pulses and tuo gamma ray pulses, none 
jf which is in phase with any of the others. For PSR-1747 and PSR-1818. 
there is one radio pulse and one gamma ray pulse but in neither case are 
the two in phase. Table 4 summarizes the situation and the relative 
shases. 


TABLE 4 

Relative Phases of Radio, Optical, X-Ray, and Gamm Ray Pulses for four 

pul sars. 


Pulsar 

I 

radio 

’base 

optical 

X-ray 

gamma ray 

very energetic 
gamma rays 

0532 (Crab) 

-20° (precursor) 






0° (main pulse) 

0° 

0° 

0°' 

variable 


143° (interpulse) 

143° 

143° 

143° 


0833 (Vela) 

0° 

100° 


60° 

• 



196° 


223° 


1747-46 

0° 



57° 


1818-04 

0° 



263° 

- 


It is unlikely that any simple model can explain such a diversity, of 
facts. The fact that radiation occurs at phases other than those which 
we would ascribe to the polar caps (i.e. the phase of the radio pulse), 
indicates that some of the radiation occurs elsewhere in the 
magnetosphere. This is also supported by the fact that high-energy 
gamma rays produced at the polar cap would be eliminated by the paii — 
production process. The Crab pulsar then appears to be anomalous in the 
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fact that pulses at all wavelengths are in phase* It is quite possible 

hat the Crab pulsar is the one pulsar in which we observe radio 

fnission from the light cylinder (or force balance) region* In that 
ase, all the radiation from the Crab is produced in one area in the 
uter magnetosphere with the possible exception of the precursor* In 
act, the precursor may be the polar-cap radio pulse that we observe In 
11 other pulsars* The fact that two pulses symmetrically arranged 
bout a center phase are observed in both the Crab and the Vela pulsars 
(for both optical and gamma ray radiation) suggests that this radiation 
is occuring in the outer magnetosphere where the opening angle of the 
open field line region is large. If this is the case, the resultant 
radiation pattern is closer to a "fan beam" than to a "pencil beam", 

thus explaining why we see both pulses in both the Vela and Crab 

pulsars. 


4*3 CONCLUSION 

To paraphrase (and invert) Voltaire's famous aphorism about God, if 
pulsars did not exist, it would not be necessary to invent them* 
Indeed, based on our current theoretical understanding of pulsars, the 
fact of their existence seems quite remarkable* Nevertheless, we can 
come to a few conclusions: 

1. If the magnetic field lines are curved (as they will be in any 
realistic model), particle acceleration (or deceleration) must 


occur . 



2 . If we assume space-charge limited f]ow> the acceleration does 
not occur near the surface of the star> but rather it occurs 
at distances of the order of a few to hundreds of star radii 
from the surface* 

3* It is likely that the star-magnetosphere system has a net 
charge^ the effect of which is to force particles emitted from 
the polar-cap region to cross magnetic field lines and 
eventually return to the star* Thus closed current loops are 
present and the system is in a (quasi) steady state. 

There is clearly much work yet to be done to establish a clear 
understanding of pulsars. Primarily we need to understand the structure 
of the magnetosphere near the light cylinder where particle inertia 
becomes important. Once we finally understand the magnetosphere, it 
will be possible to construct believable radiation emission mechanisms. 
Once the plasma processes are understood we will finally be able to make 
a strict comparison between the theoretical models and the observed 
radiation. 

It is implicit in the nature of an astrophysicist to be optimistic 
about the possibilities of understanding the distant and mysterious 
objects in the heavens. So the work on pulsar models will go on and 
better models will be produced. The task remains a difficult one and it 
is perhaps unfortunate that the LGM theory had to be dispensed with. 
Viewed in the light of what we know today it had many attractive 
features. 
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Appendix A 

MATHEMATICAL DETAILS 


1.1 GREEN FUNCTION FOR THE SECOND ORDER PERTURBATION 

The Green function for equation (3.2.64) can be easily expressed in 
:erms of the solutions to the homogeneous equation. The equation is 
Bessel's equation of order a = ^f28 and parameter A = S/TJmax- wish to 

form a Green function for the interval 0<u<u* satisfying homogeneous 
boundary conditions. For u<u' we want G(0,u')=0. For u>u' we need 

G(u*,u')=0. We therefore write G as follows: 

g(u,u^ = c ^ (A.1) 


where u^ = min(u,u') and u^ = max(u,u'). At the point u=u' there is a 
discontinuity in the first derivative given by 


dG(u,u^) 

du 


dG(u,u^) 

du 


(A. 2) 


u 


but we also know that the derivative discontinuity is given by 

CaAW(Ja,N«) (A. 3) 

where W is the Uronskian of Jo, and Na evaluated at the point u. The 
Uronskian is given by (Cf. Stegun and Abramouitz, eqn. 9.1.16) 

W(Ja,Na) = 2 /(ttAu) (A. 4) 

It immediately follows that the normalization constant 'C' is given by 

C = ir/2a (A, 5) 

Finally, the value of 'a' is determined by the boundary condition 
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5(u*>u')=0, uhich implies 


a = “J«(Xu^)/Na(^u^) 


(A, 8) 


^.2 EXACT INTEGRATION OF THE NON-^LINEAR PROBLEM 


We start from equation (3.2.7) (repeated here for convenience) 


dY 

dr 


8itej, 


|l/2 ^ ,3/2 


V' 


I \ \ p lA 

(i^) 


(A. 7) 


We noH make the substitution cosh y = y. Equation (A. 7) then can be 

written as 


/R 

y ^ = K (^) (sinh 


(A. 8) 


where 




K = — 3" 

Vv ^ 

and we may therefore integrate both sides to find 


(A. 9) 


y r 

j (sinh dy' = K j dr (A. 10) 

0 

The right hand side is trivial in integrate; the left hand side can be 
found in Dradshten and Ryzhik (1965) {equation 2.464.5 pg . 115] with the 
result given in equation (3.2.9). 
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COMPUTER LISTINGS 

The following are listings of the programs and subroutines that were 
d in the numerical solution to the non-linear one-dimensional 
ation (see Chapter 1II» section 2.3). The subroutine "ODE" and its 
support routine "DE" were written by Shampine and Gordon. All routines 
except one are written in IBM Fortran IV (level H). The one exception, 
the subroutine "INVRT" was written in the IBM 370 assembly language. 


C 

C PROGRAM TO SOLVE THE NON-LIHEAR ONE-DIMENSIONAL DIFFERENTIAL EQUATION 

C WHICH DESCRIBES THE ACCELERATION OF CHARGED PARTICLES FROM THE POLAR CAP OF A 
C 

IMPLICIT REAL*8 (A-H,0-Z) 

REAL*8 GAM(2),WORK(150) 

INTEGER IWORK(S) 

COMMON ZKAPPA,ETA,ETA2 
EXTERNAL F 
C 

C VARIABLES USED: 

C P = PERIOD 

C B = SURFACE MAGNETIC FIELD 

C RSTAR = RADIUS OF THE STAR IN CM 

C ETAMAX = MAXIMUM VALUE OF THE CO-ORDINATE ETA 

C ETAMX2 = ETAMAX**2 

C GAM(1) = RELATIVISTIC GAMMA 

C 6AM(2) = D(GAMMA)/D(XI) 

C XJSTAR = CURRENT DENSITY*! . OE- 1 2*( 1 . -ETA**2/ETAMX2 ) 

C EPSI = ACCELERATING ELECTRIC FIELD*!. OE-6 

C 

C 

C ASSUME P=! SEC 

C B=!0**12 GAUSS 

C 

RSTAR=1 ,0D+6 

ETAMX2=2.0*3. H!5927/'3. 0D+!0 
DETA=ETAMX2/!29. 

ZKAP0=2.459D!7 
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ABSERR=0.0 

RELERR=1.0D-6 


DO 10 1=1, 129,32 
ETA2=(I-1)*DETA 
ETA=DSQRT(ETA2) 

2KAPPA=2KAPO/DSQRT{1.000-.75D6*ETA2)*C 1 . 0D0-ETA2/ETAI1X2) 
XJSTAR=1 . 0D0-ETA2/ETAMX2 

GAM(1)=1.0DO+2.59D-6^XJSTAR*»(2.0DD/3.000) 

6AM(2)=-1.6D+9*(1.0D0+2.0D“1O)*XJSTAR**(2.0DO/3.0DO) 

PSIO=DSQRT (DSQRT ( 1 . ODO-1 . 0D6*ETA2) )*1 . OD-6 
WRITE(6, 110)ETA,ZKAPPA,XJSTAR 
110 FORMATC'IINTEGRATION FOR ETA = ', F 1 2 . 6, 1 0X, 'KAPPA =' 

S, 1PD20. 12,5X,'JSTAR =', F 10. 6/8X, 'PSI 12X, 'RADIUS' ,9X, 'DELTA-R' 
©,8X,'GAmA', 10X, 

«'DGAMf1A',7X,'L0G(GAMm)',7X,'E-RSl',7X,'IFLAG'/) 

PS1MAX=.25*PSI0 
PSTEP=(PSIMAX~PS10)*1.0D-9 
PSriAX = (PSinAX-PSI0)/50.D0 
CALL INVRT(PSIO,ETA,R,SN,CSN) 

DR=1.0 

GLOG=DLOG10(GArin)) 

DELTA=1.0D-4/RSTAR 

GDELT=GAn(1)“1.0D0 

EPSI=-651 .67D-6*DSQRT(1.0DO+3.0DO*CSN)/R**3/PSl0*GAn(2) 

WRITE(6, 100)PSIO,DR,DELTA,6At1(1),GDELT,6LOG,EPSI 
MRITE(9, 103)DR,GAM(1) 

HRITECIO, 103)DR,EPSI 
URITECII, 1 03) DELTA, GDELT 
103 FORriAT(1X,2020. 12) 

100 F0RMAT(1P7D15.7,5X,I5) 

P0UT=PSI0 

PSI=PSI0 

IFLAG=1 

DO 15 J=1, 1000 
POUT=PSI+PSTEP 

CALL ODE (F, 2, GAM, PSI , POUT, RELERR, ABSERR, IFLAG,WORK, IWORK) 
6LOG=DLOG10(GAM(1)) 

GDELT=GAM(1)-1.0D0 
IF(IFLAG.LT.0)GO TO 16 
CALL INVRT(PSI,ETA,R,SN,CSN) 

0R=R/RSTAR 

DELTA=(R-RSTAR)/RSTAR 

EPSI=-851 . 67D-6*DSQRT( 1 . ODO+3. 0D0*CSN)/R**3/PS1*GAM(2) 

WRITE (6, 1 00 ) PSI, DR, DELTA, GAM(I), GDELT, GL06,EPSI, I FLAG 
WRITEC9, 103)DR,GAM(1) 

WRITE! 10, 103)DR,EPSI 
URITECII, 103) DELTA, GDELT 
GO TO 18 

16 WRITEC6, 10DIFLAG, PSI, POUT 

101 F0RMAT('0***ERR0Rs I FLAG=', 14, 1 0X, 1 P2D20 . 1 2) 

STOP 

18 IFCPSMAX.LE.PSTEP)PSTEP=PSTEP+PSTEP 
IFCPSI.lt. PSI f1AX)GO TO 11 
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15 CONTINUE 


OUTPUT DATA FOR LATER PLOTTING USING 'TOP DRAWER' 

11 WRITE(9,104) 

WRITEdO, 104) 

WRITEdl, 104) 

104 F0RMAT(1X,'J0IN') 

10 CONTINUE 
STOP 
END 

SUBROUTINE F (PSI , GAM, DGAtl) 

THIS SUBROUTINE IS USED BY ODE 

IT DEFINES THE DIFFERENTIAL EQUATION TO BE INTEGRATED 

IMPLICIT REAL*8(A-H,0-Z) 

REALMS GAM(2),DGAM(2) 

COMMON ZKAPPA,ETA,ETA2 
DGAM(1)=GAM(2) 

CALL INVRT(PSI,ETA,R,SN,CSN) 

DGAM(2)=ZKAPPA*R**3*PSI**2/DSQRT( 1 . ODO- . 75D0*R*ETA2 )*GAM( 1 )/ 
EDSQRT(GAM(1)*^2-1.OD0)+6AM(2)/PSI 
RETURN 
END 


Subroutine "INVRT" is used to convert values of 7? and f to values of 
r and 8. It uses Newton's method to solve first for the value of r and 
then determines sin^(0) and cos^(8) from the definitions of t) and 


* 

* SUBROUTINE INVRT (PSI , ETA, R, SN, CSN ) 

* 

INVRT START 

USING *,15 
B START 
DC X'5',CL7'INVRT 

* 

SAV DS 18F 

* 

START EQU * 

STM 14,12, 12(13) SAV REGS 

ST 13,SAV+4 SAVE ADDRESS OM MY SAV AREA 
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SKIP1 

SKIP2 


* 


* 

ROK 


* 

^RETURN 
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LA 12,SAV 

LOAD ADDRESS OF MY SAVE AREA 

ST 12,8(13) 

STORE IN CALLING SAVE AREA 

LR 13,12 

MAKE 13 THE BASE REG. 

DROP 15 


USING SAV, 13 


LM 2, 6, 0(1) 

LOAD ARGUMENT ADRESSES 

LD 0,0(3) 

LOAD ETA 

LTDR 0,0 

TEST ETA 

BZ ZERO 

IF 0 GO TO ZERO 

MDR 0,0 

SQUARE ETA 

STD 0,ETA 

SAVE ETA**2 

LO 2,0(2) 

LOAD PSI 

STO 2,PSI 

SAVE IT 

LH 2,PSI 

LOAD UPPER HALF OF PSI 

SH 2, ETA 

SUBTRACT UPPER HALF OF ETA 

BNP SKIP1 

IF NOT POSITIVE THEN SKIP OVER DIVIDE STEP 

LD 0,=D'1.0' 

LOAD 1.0 

DDR 0,2 

1.0/PSI IN FREG 0 

B SKIP2 

SKIP AROUND THE LOAD RSAVE 

LD 0,RSAVE 

LOAD LAST VALUE OF R 

MDR 2,0 

MULTIPLY R*PSI 

LDR 4,2 

LOAD R*PSI INTO FREG 4 

MDR 4,2 

(R»PSI)**2 

MDR 4,2 

(R*PSI)**3 

LDR 6,4 

SAVE R*PSI**3 IN FREG 6 

MDR 4,2 

(R*PSI)**4 

MD 4,=D'3.0' 

3.*(PSI*R)**4 

AD 4,=D'1.0' 

3.*(R*PSI)**4+1.0 

MO 6,=D'4,0' 

4*(R*PSI)**3 

MD 6,PSI 

4*(R*PSI)**3*PSI 

AD 6, ETA 

4*(R*PSI )**3*PSI+ETA**2 

DDR 4,6 

DIVIDE TO GET NEW R 

SDR 0,4 

RSAVE-R 

LPDR 0,0 

DABS(RSAVE-R) 

STD 0, TEMPI 

SAV.E TEMPI (RSAVE-R) 

STD 4,TEMP2 

SAVE TEMP2 (R) 

LH 8,TEMP2 

LOAD UPPER PART OF TEMP2 

SH 8, TEMPI 

SUBTRACT UPPER PART OF TEMPI 

CH 8,=X'0B00' 

COMPARE WITH EXPONENT =11 

BH ROK 

IF TEMP2-TEMPD0B00 OK TO GO ON 

LDR 0,4 

LOAD NEW GUESS FOR R INTO FREG 0 

LD 2,PSI 

LOAD PSI INTO REG 2 

B SKIP2 

LOOP AGAIN 

STD 4,RSAVE 

SAVE R 

STD 4,0(4) 

RETURN R TO CALLING PROGRAM 

MD 4, ETA 

R*ETA**2 

STD 4,0(5) 

RETURN SIN**2 

SD 4,=D'1.0' 

SUBTRACT 1.0 

LPDR 4,4 

LOAD POSITIVE TO GET COS**2 

STD 4,0(6) 

RETURN CSN 


TO CALLING PROGRAM 
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RETURN L 13,SAV+4 LOAD ADDRESS OF SAVE AREA 

LM 14,12,12(13) RESTORE REGS 

MVl 12(13), X'FF' NORMAL RETURN 

BR 14 RETURN 

* 

* IF ETA=0 THEN R=1/PSI,CSN=1,SN=0 


* 


ZERO 

LD 0,=D'1.0 
DD 0,0(2) 
STD 0, RSAVE 
STD 0,0(4) 
SDR 0,0 
STD 0,0(5) 
LD 0,=D'1.0 
STD 0,0(6) 

B RETURN 

* 


* 

STORAGE 


DS 00 

RSAVE 

DC D'0.0' 

PSI 

DS D 

ETA 

DS D 

TEnPl 

OS D 

TEMP2 

DS D 


END 


DIVIDE BY PSI 

SAVE R 

RETURN R 

ZERO FREG 0 

RETURN SN=0 

LOAD 1.0 

RETURN CSN=1.0 

RETURN TO CALLING PROG 
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